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Abstract. - A resistor-network picture of transitions is appropriate for the study of energy 
absorption by weakly chaotic or weakly interacting driven systems. Such "sparse" systems reach 
a novel non-equilibrium steady state (NESS) once coupled to a bath. In the stochastic case there 
is an analogy to the physics of percolating glassy systems, and an extension of the fluctuation- 
dissipation phenomenology is proposed. In the mesoscopic case the quantum NESS might differ 
enormously from the stochastic NESS, with saturation temperature determined by the sparsity. 
A toy model where the sparsity of the system is modeled using a log-normal random ensemble is 
analyzed. 
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The study of systems with non-equilibrium steady state 
(NESS) has become active in recent years [lr[6], in- 
volving various generalizations of linear response theory 
(LRT) and of the associated fluctuation-dissipation rela- 
tion (FDR) [7}{13]. The paradigm for NESS (FigjT^i) is a 
system that is coupled to two equilibrated reservoirs, "A" 
and "B", which are characterized by different tempera- 
tures Ta and Tb- Hence the NESS of the system is not 
canonical, and it cannot be characterized by a well-defined 
equilibrium temperature. A particular case of special in- 
terest is obtained if one reservoir (call it "A") is replaced 
by a stationary driving source, while the relaxation is pro- 
vided by a bath (call it "B") that has some finite tempera- 
ture Tb- This is still the same paradigm because formally 
the driving source " A" can be regarded as a bath that has 
an infinite temperature T4 = 00. 

Stochastic modeling.— The simplest modeling of 
NESS is obtained by considering a system that has en- 
ergy levels {E n } with transition rates 



W nm = wi 



2wL 



l + c (E n -E m )/T B 



(1) 



where w e nrn and w@ m are the elements of symmetric ma- 
trices. The first term describes the transitions that 
are induced by the 7^=00 driving source. The second 
term describes the bath induced transitions, with ratio 
e (E n ~E m )/T B f n ^ m transitions, as required by detailed 
balance considerations. The dynamics of the population 
probabilities p — {p n } is described by a rate equation 



dp n 
dt 



/ Y'^nmPrn Y*mnPn\ 



(2) 



This equation can be written schematically as p — Wp, 
see remark [a]. The steady state is determined from the 
matrix equation Wp = 0. In the presence of driving the 
detailed balance is disturbed leading in general to a non- 
canonical NESS. 



Sparse systems .— In recent studies 14 -16 



our in- 
terest was focused on a class of driven systems for which 
the matrix {w e nm \ is sparse. By this we mean that the 
transition rates are characterized by a log wide (say log 
normal) distribution. In other words, the majority of el- 
ements are small, while the large elements constitute a 
small fraction, sCl. A system of current experimental 
interest is described in 
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an optical billiards with vi- 
brating walls, where the energy absorption rate (EAR) of 
the cold atoms is affected by the sparsity of the pertur- 
bation matrix, which is controlled either by the degree 
of chaoticity or by the strength of the inter-atomic inter- 
actions. Yet a simpler example concerns the absorption 



of low frequency irradiation by small metallic grains 15 



where the transitions are predominately between neigh- 
boring energy levels, and the sparsity is determined by 
the level spacing statistics. 

Beyond LRT .— In the absence of a bath the rate 
equation p = Wp generates diffusion in energy space. In 
this context it is useful to picture the levels n as "sites" 
in space, or as the "nodes" of a network, and the vf nm 
as "connectors" (Fig.ni}). The calculation of the diffu- 
sion coefficient D is exactly as the calculation of electrical 
conductivity. For example, when connecting N nodes in 
series (as in the Chain model to be discussed later), the 
"conductivity" is D = [(1/N) £ n (l/«> n )] _1 - The ada P" 



p-1 



Daniel Hurowitz and Doron Cohen 



tation of the resistor network picture to the calculation 
of D is termed semi-linear response theory (SLRT) be- 
cause D is a semi- linear function of the couplings w nm . 
This means that D[cw] = cD[w], but there is no additiv- 
ity, D[w + w'] 7^ D[w] + D[w']. Due to the sparsity, the 
result is very similar to that of percolating or disordered 
resistor-networks [17] . 

Outline.— Considering the coupling of a "sparse" sys- 
tem to a bath, our expectation is to have, as the driving 
becomes stronger, a crossover from an LRT canonical-like 
NESS to a novel non-canonical NESS, with the possibility 
of remarkable quantum-mechanical fingerprints. Specifi- 
cally our objective is to calculate the energy absorption 
rate (EAR) of a driven "sparse" system, taking into ac- 
count the presence of a surrounding environment. Be- 
low we (i) introduce the FDR phenomenology of calcu- 
lating the EAR, which requires a notion of effective NESS 
temperature; (ii) demonstrate this phenomenology for the 
simplest toy model, obtaining explicit expressions for both 
T and D; and (iii) discuss the quantum case, highlight- 
ing the existence of a saturation temperature T x that is 
determined by the sparsity. 

FDR phenomenology.— In the standard textbook 
presentation it is assumed that the system reaches a 
canonical- like state with a well-defined temperature T sys . 
The d riving i nduces diffusion with coefficient D in energy 
From E — J2 n E n p n , substituting p = Wp, 
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it follows (see appendix) that the EAR is 






E„ 



D 



with 



^ 



D [LRT] = -^w £ nm {E n -E. n 



(3) 



(4) 



where the overline indicates canonical averaging over the 
initial state. In complete analogy it is straightforward to 
show (see appendix), that the rate of cooling due to the 
interaction with the bath can be written as 



Q 



^2(E n - E m )w^ m p m 



Db 

T b 



Di 

T- 



(5) 



where the first term is due to the imbalance of upward and 
downward transition rates, while the second term is due 
to the non-uniformity of the probability distribution in en- 
ergy space. The net rate of energy increase is E = W — Q. 
At steady state W = Q, so a phenomenological determi- 
nation of the steady state temperature T aya is possible. 

The essence of the above FDR phenomenology is the 
same as in Einstein's relation: the dissipation is related to 
the diffusion (in LRT the latter is determined by the fluc- 
tuations, e.g. the velocity-velocity correlation in Einstein's 
relation, hence the terminology) . Our purpose below is to 
generalize the FDR. We emphasize in advance two issues: 



(a) The NESS might be non-canonical, so we have to de- 
fine its effective temperature; (b) The diffusion coefficient 
is not necessarily determined by LRT. 

The Chain model.— It is best to clarify the determi- 
nation of T sys and D by considering the simplest example 
of a "chain" with nearest-neighbor transitions only. With 
simplified indexing, Eq. (fijl for (n— 1) 4^> n transitions is 
written as w n + 2u> /3 /(l + e ±A °/ Ts ), where Ao is the level 
spacing, and n = 1, 2, • • •. In contrast to w° , which is the 
same for all transitions, the rates w n are characterized by 
a logarithmically wide distribution. 

In the numerical example we have N—25 levels, with 
equal level spacing Ao=l. The bath temperature is 
Tb=10. The bath induced transition rates are taken as 
w^ =0.1. The driving induced transition rates are log 
normally distributed. Namely, w n = exp(x„), where the 
x n have a Gaussian distribution with average \x and dis- 
persion a that are determined such that the driving in- 
tensity is w^ = e 2 , and the sparsity 16 is s = exp(— a 2 ). 
The value s ~ 1 means that all the elements are compara- 
ble and well represented by their average. Sparsity means 
s< 1, for which the median differs by orders of magnitude 
from the algebraic average. 

From the NESS equation Wp = we determine the oc- 
cupation probabilities p n as in FigJ2l and then calculate 
the EAR via either Eq.([f|) or Eq.([5| as in Fig[3j In the ab- 
sence of driving, the steady state is canonical with a well- 
defined temperature Tb- In the presence of driving, the 
state is generally not canonical. Consequently, we can for- 
mally associate a different microscopic temperature T nm 
for each pair of coupled levels via the defining equation 



Pu 
Pn 



exp 



E n — E n 
T 



(6) 



For the Chain model we use the simpler indexing T n , and 
the NESS equation Wp — takes the form p n /p n -i = 
Wn,n-i/Wn-i, n - Assuming A <S Tb we deduce that the 
microscopic temperature of the nth transition is given by 
the expression 



T n = 



w 



Tb 



(7) 



Effective NESS temperature.- We define the NESS 
temperature T sys such that the phenomenological FDR 
Eq.([5]) still holds. For the Chain model the bath induced 
diffusion is Db = w^Aq, and it is straightforward to show 
(see appendix) that T sys should be defined as the harmonic 
average over T„, i.e. 



T 



■«,''- 



Tn 



(8) 



For numerical results see (Fig.Hl-5). As the driving be- 
comes stronger, the temperature becomes higher, with the 
asymptotic behavior T Bys ex e 2 as implied by Eq.®. 
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LRT-SLRT crossover.— We continue with the 
stochastic Chain model and find from Eq.((3]), substituting 
Eq. (JTl) in Eq. ( 25 ) of the appendix, that the EAR is given 
by the expression 



W 



*9- 



Db 
T b 



D 



(9) 



We define the driving induced diffusion D such that the 
phenomenological FDR still holds: 



D = 



1 



;0 



w p + w T , 



n -1 



(10) 



In the limit of very weak and very strong driving we re- 
spectively get 



_D[LRT] 
Z/[SLRT1 



~A 2 



[1/W r , 



-1a2 

^0 



(11) 

(12) 



The LRT result could be obtained from the traditional 
Kubo formalism, while the SLRT prediction reflects a 
network that consists of connectors that are connected 
in series. Note that if all the w n are comparable, then 
D ss _D[lrt] sa D[slrtj. But if the w n have a log- wide dis- 
tribution, the agreement with LRT is achieved only if the 
w n are all much smaller than w^ . For strong driving, both 
D and T syB are ex [1/wn]" 1 , and hence oc e 2 , as expected 
from the SLRT resistor-network phenomenology. In this 
limit their ratio approaches the bath limited value 



W 



Db 
T b 



(13) 



which is implied by Eq.([5]). 

Quantum modeling.— It should be clear that SLRT 
applies whenever the transport is modeled using a resistor 
network. Thus it might have applications, e.g., in statis- 
tical mechanics and biophysics. But the original motiva- 
tion for SLRT came from mesoscopics, where the quantum 
nature of reality cannot be ignored. In this context the 
Hamiltonian contains a driving term f(t)V, and the tran- 
sition rates w e nm oc \V nm \ 2 between levels are determined 
by the Fermi-Golden- Rule (FGR). For weakly chaotic or 
weakly interacting systems V nm is typically sparse and tex- 
tured. Assuming that /(£) has zero average, and correla- 
tion function (f(t)f(t')) = e 2 S(t — if), the rate equation 
Eq.S is replaced by the quantum Master equation 



dp .r^ , £ 2 

Tt = - l[n °' p] -J 



V,{V,p]} + W p p 



(14) 



where the second and third terms correspond to the driv- 
ing source and to the bath as in Eq.ll]). This equation is 
of the form p = Wp. It is identical with p = Wp of Eq.pl) 
if the off-diagonal terms are ignored, where 



« ! „ 



e 2 \V n 



Note that W induces dephasing of the off-diagonal ele- 
ments, which we assume to be minimal in the "quantum" 
case: see technical details below. 



Some technical details with regard to Eq. ( 14 1 are in 



order (can be skipped in first reading). In the energy basis, 
"Ho is a diagonal matrix with energy levels E n . The state 
of the system is represented by the probability matrix, 
which can be rewritten as a column vector p 1— > (jp n ; /9„ M ), 
composed of the diagonal probabilities and the off-diagonal 
coherences. Consequently, the master equation takes the 
form p = Wp, with the super operator 



W 



W At 

A W 1 - 



(16) 



?\V nm \\ 
The definition of A is implied by the second term in 



The matrix W is given by Eq. (hi) with w z nm = e 2 \V nm \ 



Eq.(14). In particular we note that 



w. 






= s 2 V nv V„ 



y fini 



IP (17) 

for i/,nj^n (18) 



We use the common notations A„ M = E v —E^, and 7^ = 
(e 2 /2)[(V 2 ) uu +(V 2 ) tipi ]. For simplicity, we assume that 
the bath induced dephasing 7,3 = ufi + "f v is the same for 
all the coherences. 



In the absence of driving Eq.( 14 ) is the well known Pauli 



master equation leading to canonical equilibrium. In the 
numerics we assume in the "quantum" case minimal de- 
phasing, which means j v — 0. More generally we may 
have an extra "pure dephasing" effect which is represented 
by If > as m the familiar NMR context. 

The stochastic (FGR) picture applies if the coherences 
become negligible. By inspection, it might happen for two 
reasons. One possibility is that there is strong pure de- 
phasing effect (7 V » w £ ) that suppresses the coherences. 
Then one simply recovers the stochastic rate equation for 
the occupation probabilities. The second possibility is to 
have w e much smaller compared with the level spacing Ao. 
The latter is the traditional assumption in atomic physics, 
and can be regarded as "microscopic circumstances" . But 
in "mesoscopic circumstances" Ao might be small, and the 
validity of FGR is not guaranteed. 

Quantum NESS.- The NESS is obtained by solving 
the equation Wp = 0. Looking at Fig.[4]-5 we see that 
in the "quantum" case there is a saturation temperature 
Tb < Too < 00 that depends on the sparsity s. Conse- 
quently there is a premature saturation of the EAR at 



W^ = 



Db 

T B 



Db 

T 



(19) 



as seen in Fig[3] Note (El . Only in the non-sparse limit do 
we recover quantum-to-classical correspondence. 

One would like to understand how X^ emerges. For this 
purpose we regard Eq. ( 14 1 for p nn > as a Fokker-Planck 
equation, where n is the momentum. In the absence of 
sparsity, V nn < can be interpreted as the matrix represen- 
(15) tation of the position coordinate, and its eigenstates are 
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extended in n. But if s <C 1, then V nn ' is like off-diagonal 
disorder, and its eigenstates r become localized in n, with 
energies (E) r that are no longer identical. If the driving 
is very strong, the master equation implies that the NESS 
becomes a mixture of the eigenstates r with weights 



p r ex exp 



(E)r 



TL 



(20) 



In the sparse limit the n and the r bases essentially coin- 
cide, and accordingly T^ ~ T mix ~ Tg. This is confirmed 
numerically in FigjHJ 

In order to further illuminate the variation of the sat- 
uration temperature as a function of s = exp(— a 2 ), we 
have added in Fig[6]a lower bound estimate for X^ which 
is deduced as follows [cj: Considering an energy range 
A(E n ) of interest, one observes that as a is decreased, the 
(E) r occupy a smaller range A(E r ), and eventually in the 
non-sparse limit they are quasi-degenerate. The p r distri- 
bution is characterized by some T m i X > Tg . Consequently 
the stretched distribution p n is characterized by a higher 
temperature X^ > [A(E n ) / A(E r )]Ts ■ Hence the satura- 
tion temperature increases for smaller a, and diverges in 
the non-sparse limit. 

Experimental Signature.— The thermodynamic defi- 
nition of temperature which we have adopted via Eq.([5]) is 
related to the heat flow between two bodies: it is the same 
concept which is reflected in the phrasing of "the zeroth 
law of thermodynamics" and in the associated definition 
of empirical temperature. It should be contrasted with 
an optional statistical mechanics perspective that relates 
the temperature to the fluctuations, as in the paradig- 
matic Brownian motion studies of Einstein. Once T Bya of 
Eq. (pi) is determined through a heat flow measurement, 
it is possible to determine experimentally also the D in 
Eq.© through an EAR measurement. Then it is possible 
to explore the LRT to SLRT crossover of this absorption 
coefficient Eq. pT])-Eq. (12). This crossover is reflected in 
having a wide range (many decades) over which the EAR 
varies as in Fig|3] It is important to realize that the un- 
limited increase of temperature in the classical case does 
not have by itself a reflection in the EAR, because clas- 
sically the EAR always saturates to a bath limited value 
as determined by Eq. (13). In contrast to that, quantum 



mechanically there is a premature saturation, as expected 
from Eq. (19), reflecting the sparsity of the system. This 



effect might be useful for the detection of a classical to 
quantum transition with non-equilibrium measurements. 
Discussion.— The "sparse" NESS resembles that of a 
glassy system 11 : In both cases, one has to distinguish 
between "microscopic" and "macroscopic" time scales, 
rates, and temperatures. Eq. ([3| and Eq. ^ establish 
a diffusion-dissipation relation involving a "macroscopic" 
temperature T that might be much lower compared with 
the microscopic temperatures T n . The diffusion is driven 
by the fluctuations of the sources, but it is not the LRT 
(Kubo) formula which should be used in order to deter- 
mine D. but rather a resistor network SLRT calculation. 



A closely related chain of works regarding NESS, con- 
cerns mixed phase-space of periodically driven systems pl- 
6 , where the problem is reduced to the study of a stochas- 
tic rate equation. Our work differs in three respects: 
(1) The sparsity may arise even for quantized chaotic non- 
mixed systems, implying a glassy type of NESS; (2) We 
assume a stationary driving source instead of a strictly 
periodic driving; (3) The master equation approach has 
allowed us to consider novel mesoscopic circumstances in 
which the quantum NESS differs enormously from the 
stochastic prediction. 

The influence of quantum coherence on the NESS is 
remarkable. Due to the localization of the eigenstates in 
energy space, we found that for strong driving the temper- 
ature saturates to a finite value that reflects the sparsity 
of system. This should be contrasted with the traditional 
prediction of unbounded temperature. 

Appendix.— For the convenience of the reader we con- 
centrate here the technical details that concern the deriva- 
tion of the FDR phenomenology. The energy of the system 
is E = J2 n P™ E n, and its rate of change is E = J2 n E nt>n- 
The equation for p n includes a driving source term and a 
bath term. Accordingly we write E = W — Q, where the 
two terms are interpreted as the rate of heating due to the 
driving (equals the EAR), and the rate of cooling due to 
the bath (in steady state equals the EAR too). From the 
master equation it follows that the expression for Q, both 
in the stochastic and in the quantum case, is given by the 
first equality of Eq. |5]) . This expression can be written as 
the sum of a term that originates from the asymmetry of 
w@ m , and a term that originates from the non-uniformity 



of the p n . Defining p r , 



w p -w 13 = 

nrn ^mn 



(p n +p m )/2 we get 



Pn Pra 



2 tanh - 
2 tanh ( - 



2T B 

£j - L nrn. 



Pn 



(21) 
(22) 



At high temperatures one can approximate the tanh() by 
linear functions leading to 



Q 



1 



if. 



r> / , Pnm rri K-^n -"r? 



1 v^ - < m 

/ ,P1 



(E n —E„ 



(23) 
(24) 



The first expression is identified as Db/Tb, where Db is 
the bath induced diffusion coefficient. The second expres- 
sion is used to define the effective system temperature T Bys , 
such that it takes the form —Db/T byb . 

In the stochastic case the effect of the driving can be 
treated using the same procedure. The expression for W 
is given by the first equality of Eq. pi , that is analogous 
to Eq.([5j). Taking into account that w e nm unlike w^ m is a 
symmetric matrix one obtains, in analogy to Eq.(24), 






-(En—Em) 



(25) 
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If the state is strictly canonical with all the T nm equal 
the same number T sya , then W = D/T sys where D is given 
by Eq. Q . This is the standard LRT expression for the 



diffusion coefficient, leading to the linear result Eq.( 11 ) in 
the case of near-neighbor transitions. 

More generally T sys is the effective temperature, and the 
equation W = D/T aya is used to define the effective diffu- 



sion coefficient D, leading in the stochastic case to Eq.( 10 ) 



This agrees with the linear result Eq.( 11 ) for weak driving 



and with the semi-linear result Eq.(12| for strong driving. 
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Bath Driving 



Fig. 1: (a) Block diagram of the model system. Both the driv- 
ing source (A) and the bath (B) induce fluctuations that drive 
transitions. These fluctuations are characterized by a power 
spectrum S(u>). The driving source can be formally regarded 
as a bath that has an infinite temperature Ta — oo, such that 
§a(— w) = Sa(u). In contrast to that Tb is finite, and ac- 
cordingly Sb{—u)/§b{u) = exp(w/TB). (b) Illustration of 
the Chain model. In the absence of a bath the driving induces 
transitions (red arrows) between levels E n of a closed system, 
leading to diffusion in energy space and hence heating. The 
diffusion coefficient D can be calculated using a resistor net- 
work analogy. A NESS is reached due to the presence of a heat 
bath (blue arrows) that favors downward transitions. 
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Fig. 2: The NESS occupation probabilities p n are plotted vs 
E n in the stochastic (blue) and quantum (red) cases. In the 
latter (quantum) case we plot also the occupation probabilities 
p r of the V eigenstates versus {E) r - The sparsity is s = 10 , 
and e = 9.3. We observe that the effective temperature pre- 
dicted by the quantum master equation is lower compared with 
the stochastic approximation. 
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Fig. 3: The EAR W for the stochastic (solid blue) and for the 
quantum (dashed red) NESS. The sparsity is s — 10 -5 . The 
vertical lines are plotted at values of e for which the stochastic 
picture predicts a crossover: i.e. the e values for which w^ and 
[1/iUn]" 1 equal w". 



Fig. 5: The effective NESS temperature T 5ys in the stochastic 
case (upper panel) and in the quantum case (lower panel) is 
imaged for additional values of the sparsity. (Color online) Blue 
represents the bath temperature T = 10, while red corresponds 
to T = 50. For quantitative dependence see FigJ4] 
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Fig. 4: The effective NESS temperature T sys versus the driv- 
ing intensity. Solid (blue) lines are for the stochastic NESS, 
while dashed (red) lines are for the quantum NESS. The dot- 
ted (green) line represents the temperature Tb of the bath. 



Fig. 6: The dependence of Too on the width a of the log-normal 
distribution. Note that the sparsity is s = exp(— a 2 ). We con- 
firm that Too is bounded from below by [A(E n ) / A(E r )]TB 
(dashed red line), and tends to Tb in the sparse limit. Here 
A(E n ) = 25 is the width of energy window in this numerical 
test, while A(j5" r ) is the length of the interval that contains the 
energies (E) r . 
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